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The asymptotic probability density function of nonlinear phase noise, often 
called the Gordo n-MoUenauer effect, is derived analytically when the number 
of fiber spans is very large. The nonlinear phase noise is the summation 
of infinitely many independently distributed noncentral chi-square random 
variables with two degrees of freedom. The mean and standard deviation of 
those random variables are both proportional to the square of the reciprocal 
of all odd natural numbers. The nonlinear phase noise can also be accurately 
modeled as the summation of a noncentral chi-square random variable with 
two degrees of freedom and a Gaussian random variable. @ 2008 Optical 
Society of America 
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When optical amplifiers are used to compensate for fiber loss, the interaction 
of amplifier noise and the Kerr effect causes phase noise, often called the Gordon- 
Mollenauer effect or nonlinear phase noise Nonlinear phase noise degrades both 
phase-shifted keying (PSK) and differential phase-shift keying (DPSK) systemJ^ 
that have renewed attention recently.'^'^ Usually, the performance of the system is 



estimated based on the variance of the nonhnear phase noise .^However, the nonhnear 
phase noise is not Gaussian nois^ and the variance is not sufficient to characterize 
the system. The probabihty density function (p.d.f.) is required to better understand 
the system and evaluates the system performance. This letter provides an analytical 
expression of the asymptotic p.d.f. for the nonlinear phase noise when the amplifier 
noise is modeled as a distributed process for a large number of fiber spans. The 
characteristic functions are ffist derived analytically as a simple expression and the 
p.d.f is the inverse Fourier transform of the corresponding characteristic function. The 
asymptotic p.d.f. can be accurately applied to system having more than 32 spans. 
For an A^-span fiber system, the overall nonlinear phase noise i^ 

0LN = 7^cfr {|A + nip + |A + m + nsp + ■ ■ ■ + \A + m + ■ ■ ■ + n^vp} , (1) 

where A is a real number representing the amplitude of the transmitted signal, Uk, k = 
1, . . . ,N, are independent identically distributed (i.i.d.) complex zero-mean circular 
Gaussian random variables as the optical amplifier noise introduced into the system 
at the kth fiber span, jLes is the product of fiber nonlinear coefficient and effective 
fiber length per span. 

With large number of fiber spans, the summation of ^ can be replaced by 
integration as 

^LN = f^ r \A + Siz)\Mz, (2) 
Jo 

where S{z) is a zero-mean complex value Wiener process or Brownian motion of 
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E {S (zi) S* {z2)} = (T^ mill (2:1, 2:2) and k, = N'-fLcs/L is the average nonlinear coeffi- 
cient per unit length. The variance of cr^ = Na\^-^/L is the noise variance per unit 
length where = a\^^, k = 1, . . . , is noise variance per amplifier. 

The p.d.f. is derived for the following normalized nonlinear phase noise 

0= f'\p + bit)\'dt, (3) 
Jo 

where b{t) is a complex Wiener process with an autocorrelation function of 

Rb{t, s) = E{b{s)b*{t)} = min(t, s). (4) 

Comparing the integrations of ^ and Q , the normalized phase noise of Q is scaled 
hj (f) = Lag(l)i^t^/ K, t = z/L is the normalized distance, b{t) = S{tL) /ag/^/L is the 
normalized amplifier noise, and p = A/as/VL is the normalized amplitude. The 
optical signal-to-noise ratio (SNR) is = ^^/(Laf) = /{Naj^^^). 

The Wiener process of b{t) can be expanded using the standard Karhunen-Loeve 
expansion of 7, §6-4] 

00 

b(t) = (JkXk^pkit), (5) 

k=l 

where Xk are i.i.d. complex circular Gaussian random variable with zero mean and 
unity variance, al,iljkit),0 < t < 1 are the eigenvalues and eigenfunctions, respec- 
tively, of the following integral equation, 

Mi) = (^l f' Rb{t,s)Ms)ds, (6) 
Jo 



3 



with boundary condition of ipk{0) = 0. The eigenfunctions of ipk{t) are orthonormal 



1 1 k = l 

Mt)i'imt={ ■ (7) 

k=^l 



Substitute the correlation function of (jH) into the integral equation of (jHI), we get 

i'k{t) = crl [ sipkis)ds + alt [ i'k{s)ds. (8) 
Jo Jt 

Take the second derivative of both sides of (jH)) with respect to t, we get 



^iMt) (9) 



with solution of ilj{t) = \/2sm{t/ak)- Substitute into (0) or (jH)), we find that 



{2k - 1)71^ 



(10) 



Previous studieJ^ are equivalent to the Karhunen-Loeve transform of finite num- 
ber of random variables of based on numerical calculation. While the eigenvalues 
of the covariance matrix corresponds approximately to al of (fTUj) . the eigenvectors 
always require numerical calculationsP The assumption of a distributed process of 
(121) can derive both eigenvalues and eigenfunctions of (fTIHl analytically. 

Substitute (0) with (fTUI) into the normalized phase of ©, because Jq sm{t/ak)dt = 
(jfc, we get 



/ + 2V2Y.al^{xk) + E '^ll^kl'. (11) 

k=l k=l 
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where 3?(-) denotes the real part of a complex number. Because J2k'=i ^1 = 1/2 (see 
§0.234]), we get 

oo 

^ = Y.al\V2p + Xk\'. (12) 

k=l 

The random variable \ V2p + is a noncentral x-SQuare random variable with 
two degrees of freedom with a noncentrality parameter of 2p^ and a variance parameter 
of 1/2 JOl p. 44]. The normalized nonlinear phase noise is the summation of infinitely 
many i.i.d. noncentral x-square random variables with two degrees of freedom with 
noncentrality parameters of 2cr^p^ and variance parameter of cr^/2. The mean and 
standard deviation of the random variables are both proportional to the square of the 
reciprocal of all odd natural numbers. 

The characteristic function of \\/2p + x^p is (TUl p.44] 



'1, (13) 

and with mean and variance of 2p^ + 1 and 4p^ + 1, respectively. The characteristic 
function of the normalized phase of is 



(14) 



k=i^-3^^k V-J^(^kJ 
Using the expressions of jHl §1.431, §1.421], the characteristic function of (jl4j) can be 

simplified to 



(15) 



^<t>{ji^) = sec{Jjiy)exp 
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The first eigenvalue of pOj) is mucli larger than other eigenvalues. The normalized 
phase of (|TT|l is dominated by the noncentral x-square random variable corresponding 
to the first eigenvalue because of 



^ 



and 



1/2 - (2/7r)^ 



4.27, 



(16) 



cy\ + at + 



68.12. 



17) 



1/6 - (2/7r)4 

The relationship of Z^fcli c"^ = 1/6 is based on ^9, §0.234]. 

Beside the noncentral X"Square random variable corresponding to the largest 
eigenvalue of ai, the other x-square random variables of |v^p + Xfcp, A; > 1, have more 
or less than same variance. From the central limit theorem [3 §5-4], the summation 
of many random variables with more or less the same variance approaches a Gaussian 
random variable. The characteristic function of (|T^ can be accurately approximated 
by 



1 / SjzV/vr^ 

exp — — 

1 - Ajv/Ti^ \1 ^ Ajv/'K^ 

X exp jz/(2p^ 



•1 4 
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(i- 





as a summation of a noncentral x-square random variable with two degrees of freedom 
and a Gaussian random variable. While the characteristic function of (fT3j) is a simpler 
expression than that of (fTH|l . the physical meaning of (fTK|) is more obvious. 
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The p.d.f. of the normahzed phase noise of Q can be calculated by taking the 
inverse Fourier transform of either the exact (|T5|l or the approximated ()18|1 charac- 
teristic functions. Fig. ^ shows the p.d.f. of the normalized nonlinear phase noise for 
three different optical SNR of = 11, 18, and 25, corresponding to about an error 
probability of 10^^, 10~^, and 10^^^, respectively, when amplifier noise is the only 
impairment. Fig. [T] shows that the p.d.f. using the exact (fT3j) or the approximated 
(fTHj) characteristic function, and the Gaussian approximation with mean and variance 
of = + 1/2 and cr| = (4p^ + l)/6. The exact and approximated p.d.f. overlap 
and cannot be distinguished with each other. 

Fig. 121 shows the cumulative tail probabilities as a function of Q-factor. The Q- 
factor is defined a.s Q = {(f) — m^)/a^ and gives an error probability or tail probability 
of |erfc((5/"\/2) for Gaussian distribution, where erfc(-) is the complementary error 
function. Fig. |21is plotted for the case of = 18. From Fig. |2l the p.d.f. calculated 
from the exact (fT3|) or approximated ()18p characteristic function has no difference. 
The Gaussian approximation underestimates the cumulative tail probability for Q > 1 
but overestimates the cumulative tail probability for Q < 1. 

The p.d.f. for finite number of fiber spans was derived base on the orthogonal- 
ization of (Q) by iV i.i.d. random variables.!^ Fig. 01 shows a comparison of the p.d.f. 
for = 4, 8, 16, 32, and 64 of fiber spanJ^ with the distributed case of ()15p . Using an 
optical SNR of = 18, Fig. El is plotted in logarithmic scale to show the difference 
in the tail. Fig. IHl also provides an inset in linear scale of the same p.d.f. to show the 
difference around the mean. The asymptotic p.d.f. of p5|) with distributed noise has 
the smallest spread in the tail as compared with those p.d.f.'s with N discrete noise 
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sources. The asymptotic p.d.f. is very accurate for N > 32 fiber spans. 

In summary, this letter derives the asymptotic p.d.f. of nonhnear phase noise 
when the number of fiber spans is very large. Gaussian approximation based solely 
on the variance cannot use to predict the performance of the system accurately. The 
nonlinear phase noise can be modeled accurately as the summation of a noncentral 
square random variable with two degrees of freedom and a Gaussian random variable. 
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List of Figure Captions 

Fig. ^ The p.d.f. of the normahzed nonhnear phase noise for optical SNR of = 
11,18, and 25. 

Fig. 121 The cumulative tail probability as a function of Q-factor. 

Fig. El The asymptotic p.d.f. of as compared with the p.d.f. of = 4, 8, 16, 32, and 

64 fiber spans. The p.d.f. in linear scale is shown in the inset. 
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Normalized nonlinear phase noise, (j) 

Fig. 1. The p.d.f. of the normahzed nonhnear phase noise for optical SNR 
of p2 = 11,18, and 25. 
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Normalized nonlinear phase noise, 



Fig. 3. The asymptotic p.d.f. of as compared with the p.d.f. of N — 
4, 8, 16, 32, and 64 fiber spans. The p.d.f. in hnear scale is shown in the in- 
set. 
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